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ABSTRACT 

We present the PLANET photometric data set for OGLE-1998-BUL-14, a high magnification (A max ~ 
16) event alerted by the OGLE collaboration toward the Galactic bulge in 1998. The PLANET data 
set consists a total of 461 /—band and 139 V— band points, the majority of which was taken over a 
three month period. The median sampling interval during this period is about 1 hour, and the la 
scatter over the peak of the event is 1.5%. The excellent data quality and high maximum magnification 
of this event make it a prime candidate to search for the short duration, low amplitude perturbations 
that are signatures of a planetary companion orbiting the primary lens. The observed light curve for 
OGLE-1998-BUL-14 is consistent with a single lens (no companion) within photometric uncertainties. 
We calculate the detection efficiency of the light curve to lensing companions as a function of the mass 
ratio and angular separation of the two components. We find that companions of mass ratio > 0.01 
are ruled out at the 95% confidence level for projected separations between 0.4 — 2Ate, where te is the 
Einstein ring radius of the primary lens. Assuming that the primary is a G-dwarf with r~E ~ 3 AU our 
detection efficiency for this event is ~ 60% for a companion with the mass and separation of Jupiter and 
~ 5% for a companion with the mass and separation of Saturn. Our efficiencies for planets like those 
around Upsilon And and 14 Her are > 75%. 

Subject headings: gravitational lensing, dark matter, planetary systems 
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1. INTRODUCTION 

Mao & Paczyhski (1991) first suggested that planets 
could be detected in microlensing events; Gould & Loeb 
(1992) pointed out that if all stars had Jupiter-mass plan- 
ets with separations near 5 AU, then ~ 20% of the mi- 
crolensing events should exhibit detectable planetary de- 
viations, provided that events were monitored frequently 
and with moderately high precision. Current microlens- 
ing discovery teams do not generally sample frequently or 
precisely enough to detect the short-lived perturbations 
caused by planets. However, since these collaborations re- 
duce their data in real time, they are able to issue 'alerts,' 
notification of ongoing microlensing events detected before 
the peak magnification. Prompted by this alert capabil- 
ity, several other groups have formed to monitor alerted 
microlensing events more closely (GMAN, Alcock et al. 
1997; PLANET, Albrow et al. 1998; MPS, Rhie et al. 
1999a). Since only a handful of alerted events are in 
progress at any given time, they can be monitored with 
the fine temporal sampling and photometric precision re- 
quired to discover planetary perturbations. In particular, 
the PLANET (Probing Lensing Anomalies NETwork) col- 



laboration has access to four telescopes that are roughly 
equally-spaced in longitude, and thus can monitor mi- 
crolensing events almost continuously, weather permitting. 

Over the last five years, PLANET has monitored over 
100 events with varying degrees of sampling and photo- 
metric precision. Here we present photometry and analy- 
sis of one such event, OGLE-1998-BUL-14, the 14th event 
alerted by the OGLE collaboration toward the Galac- 
tic bulge in 1998. The total PLANET data set for this 
event consists of 600 data points, the majority of which 
was taken during a 3 month period starting 1 May 1998. 
The median sampling interval during this interval is about 
1 hour, with no gaps greater than 4 days. The photometric 
precision near the peak of the event, where the sensitivity 
to planets is highest, is 1.5%. These characteristics, com- 
bined with the high maximum magnification of OGLE- 
1998-BUL-14 make our data set highly sensitive to plane- 
tary perturbations. The PLANET data for the event are 
consistent with a generic point-source point-lens (PSPL) 
model. The excellent photometry and dense sampling also 
allow us to place stringent constraints on possible stellar 
and planetary companions. We also place limits on paral- 
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lax effects arising from the motion of the Earth, deviations 
arising from the finite size of the source, and the amount 
of blended light from the lens itself. These limits are then 
translated into limits on the mass of the lens. We find that, 
despite the excellent coverage and photometry of OGLE- 
1998-BUL-14, the limits on the mass of the lens are very 
weak. This indicates that it will in general be quite diffi- 
cult to obtain interesting constraints on the masses of the 
lenses giving rise to microlensing events from photometric 
data alone. 

Our study is similar to that done by the MPS and MOA 
collabora tions on the micr olensing event MACHO-1998- 
BLG-35 ( Rhie et al. 1999b ), which was a higher maximum 
magnification event (A max ~ 75) than OGLE-1998-BUL- 
14. We compare the limits on companions for OGLE-1998- 
BUL-14 to those for MACHO-1998-BLG-35 directly in § 
5.2. 

A brief introduction to the theory of microlensing is 
given in § 2. The observations and data are presented in § 
3. In § 4, we discuss known systematic effects in crowded- 
field photometry and our method of correcting for them, 
and fit the data to a PSPL model. In § 5, we search for 
the kinds of deviations from the PSPL model that would 
arise from companions to the primary lens. Finding none, 
we calculate the detection efficiency of OGLE-1998-BUL- 
14 light curve to companions as a function of the mass 
ratio and angular separation of the companion, and use 
this efficiency to place limits on possible companions to 
the primary lens. In § 6, we use several considerations to 
constrain the mass of the primary lens. In § 7, we convert 
from mass ratio and angular separation to mass and physi- 
cal separation of the companion using an assumption of the 
mass and distance to the primary lens, and compare the 
OGLE-1998-BUL-14 detection efficiencies to other meth- 
ods of detecting extrasolar planets. We conclude in § 8. 

2. BASIC MICROLENSING 
The flux of a microlensing event is given by 



F(t) = F s A{t) + F B , 



(1) 



where Fg is the unlensed flux of the source, A(t) is the 
magnification as a function of time, and Fb is the flux of 
unresolved stars not being lensed, which may include light 
from the lens itself. For a PSPL model, the magnification 
is 



AW)\ = 



(2) 



where u(i) is the angular separation of the source and the 
lens in units of the angular Einstein radius 9 E defined by 
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where M is the mass of the lens, and Dls 5 Ds> -Dl are the 
lens-source, observer-source, and observer-lens distances, 
respectively. This corresponds to a physical distance at 
the lens plane of 



r E = E D h ~ 3 AU 
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For the scaling relation on the far right of equations (Q) 
and (||), we have assumed D$ = 8 kpc and D E = 6.5 kpc, 
typical distances to the lens and source for microlensing 
events toward the bulge. For rectilinear motion, 
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where to is the time of maximum magnification, u Q is the 
minimum impact parameter of the event in units of 9 E , and 
tE is the Einstein time scale, a characteristic time scale of 
the event defined by 



t E = 
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Here v is the transverse velocity of the lens relative to the 
observer-source line-of-sight. For the scaling relation on 
the far right of equation (g) we have assumed a transverse 
velocity of v = 130 kms -1 , and we have again assumed 
Dl = 6.5 kpc and D$ — 8 kpc. 

A PSPL fit to an observed data set is a function of 
3 + 2 Ni parameters: tE, uq, to, and one source flux Fg 
and one blend flux Fb for each of Ni light curves taken in 
different sites or in different bands. 

3. OBSERVATIONS 

PLANET observations of OGLE-1998-BUL-14 were 
taken in two broad-band filters at five sites using seven 
different detectors. The five sites are the CTIO 0.9m 
and the Yale-CTIO lm in Chile, the SAAO lm in South 
Africa, the Canopus lm in Tasmania, and the Perth 0.6m 
near Perth, Australia. Canopus data prior to HJD' = 
HJD - 2450000.0 = 975.0 were taken with a different de- 
tector than those taken afterwards; we will refer to these 
data sets as Canopus A (HJD' < 975.0) and Canopus B 
(HJD' > 975.0), respectively. SAAO data were taken in 
three segments: during the period from HJD' = 976.0 to 
HJD' = 980.0, a different detector was used than prior 
to HJD' = 976.0 or after HJD' = 980.0, when the orig- 
inal detector was reinstalled on the telescope. Since dif- 
ferent detectors have different characteristics that can af- 
fect the photometry, we will treat these as independent 
light curves. In addition, because the SAAO data prior 
to HJD' = 976.0 and after HJD' = 980.0 are substantially 
offset photometrically despite being taken with the same 
telescope, detector and filters, these are also treated as 
independent light curves. We refer to these as SAAO A 
(HJD' > 976.0), SAAO B (976.0 < HJD' < 980.0) and 
SAAO C (HJD' > 980.0). All independent light curves 
have both / and V band photometry, except for the CTIO 
0.9m and Canopus B data, which do not have V^-band 
photometry. 

The entire data set consists of 461 /-band and 139 V- 
band data points forming a total of 14 independent light 
curves. The number of data points per light curve is 
given in Table 1. The data were reduced using DoPHOT 



(Schechter, Mateo, & Saha 1993), and reference stars were 
chosen to optimize photometry at each individual site. For 
details concerning the reduction, see Albrow et al. (1998). 

The PLANET data set for OGLE-1998-BUL-14 is one 
of the best sampled light curves to date. Over 95% of 
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Fig. 1. — Histogram of the distribution of sampling intervals (time between successive measurements) in hours for the OGLE-1998-BUL-14 
dataset. The median sampling interval is about 1 hour. 



the measurements were taken during a 3.3 month time pe- 
riod from HJD' = 945.0 to HJD' = 1045.0, corresponding 
to times — 0.2£e before the peak until 2£e after the peak. 
In Figure 1, we show the distribution of sampling inter- 
vals (time between successive measurements) in hours for 
the cleaned PLANET OGLE-1998-BUL-14 dataset. The 
median sampling interval is about 1 hour, or 10~ 3 of the 
Einstein ring crossing time. Furthermore, there are very 
few gaps greater than 1 day. 

Our primary results are based solely on PLANET 
data. We use the publically-available OGLE data set 1 for 
OGLE-1998-BUL-14 only to test our PSPL model and to 
derive a parallax-based constraint on the lens mass. The 
OGLE data set consists of 159 data points taken in the 
standard I filter with the 1.3 Warsaw Telescope in Chile: 
125 baseline points were taken prior to HJD' = 800 when 
the source was not being lensed, and the remaining 34 
points taken during the course of the event. For more in- 
formation on the OGLE project and the Early Warning 
System, which is used to alert microlensing events toward 
the Galactic bulge, see Udalski, Kubiak, & Szymanski 
(1997), and Udalski et al. (1994). 

1 http://www.astrouw. edu.pl/~ftp/ogle/ogle2/ews/ews. html 



4. CORRECTING FOR KNOWN SYSTEMATIC EFFECTS 

Along with the usual reduction procedures, we take ad- 
ditional steps to optimize the data quality before analyz- 
ing the light curve of OGLE-1998-BUL-14. Since plane- 
tary perturbations are expected to often have small am- 
plitudes, it is essential that any low-level systematic de- 
viations caused by observational effects be minimized in 
order to avoid spurious detections. Such effects can be 
quite common in crowded-field photometry. We describe 
in turn several systematic effects and our procedures to 
remove them. 

Light curves of constant stars in our fields often display 
residuals from the mean value for the star that are cor- 
related with seeing (or more specifically, image quality as 
measured by the FWHM of point sources), and occasion- 
ally with the sky background as well. These correlations 
seem to be a generic feature of crowded-field DoPHOT 
photometry, and presumably are present at some level in 
all light curves. We find that the magnitude (and sign) 
of the correlation depends on the site and detector, being 
quite strong for some data sets and almost completely ab- 
sent in others. Typically, these correlations with the seeing 
and background are linear in flux and hence approximately 
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Fig. 2.— The combined PLANET (open circles) and OGLE (crosses) light curves for OGLE-1998-BUL-14. In order to show all data 
on the same figure, we have plotted the total magnification for each light curve, which is given by A = (F — F&)/Fs, where F is the ob- 
served flux, and Fg and Fg are the source and blended fluxes derived by fitting the data to a PSPL model, respectively. All data prior to 
HJD' = HJD - 2450000 ~ 940 are from OGLE; the majority of the data after HJD' = 940 are from PLANET. 



linear in magnitudes, and are below 10%. Such correla- 
tions innate the overall scatter in typical light curves by 
about a factor of two, diminishing the recognizability of 
subtle deviations. The correlation with FWHM can be 
especially dangerous: coherent deviations of order a few 
percent are seen on nightly time scales, partly due to the 
correlation between seeing and airmass. Such deviations 
can easily mimic low amplitude perturbations caused by 
small-mass planetary companions. 

Based on studies of constant stars in crowded fields, 
Albrow et al. (1998) found that formal DoPHOT errors 
typically underestimate the true photometric uncertain- 
ties by a factor of 1.5 — 2. Using the formal DoPHOT 
errors for our analysis would therefore overestimate the 
significance of any anomaly being studied. Furthermore, 
observed error distributions are not Gaussian, with long 
tails toward larger values, primarily due to the seeing and 
background systematics described above. These distribu- 
tions are poorly represented by the formal DoPHOT er- 
rors. Thus, when not corrected for systematics, many light 
curves have more large (> 3cr) outliers than would be ex- 
pected from a Gaussian distribution. While it may be 
tempting simply to eliminate these outliers from the data 



set, such an approach could be dangerous since an isolated 
outlier could, in principle, be due to a short-duration de- 
viation caused by a planetary companion. Unfortunately, 
in many cases the cause of the outliers is not known and 
nearly-simultaneous photometry that in principle could be 
used to discriminate real deviations from poor-quality data 
is not available. 

Our approach is to apply a correction to the entire 
dataset prior to the analysis, as we now describe. We 
fit the entire data set to a preliminary PSPL model, in- 
cluding seeing and background correlation terms, so that 
for each observed light curve our model takes the form, 

niprccu = topspl(^) + r\Q% + Ch- (7) 

Here m prc d,i is the predicted flux in magnitudes, mpspL 
is the flux in magnitudes due to the PSPL magnification 
(Eq. [|) at time U, r\ is the slope of the seeing correlation, 
6i is the DoPHOT-reported full-width half-maximum of 
the PSF of the ith data point, £ is the slope of the correla- 
tion with background flux, and 6, is the DoPHOT-reported 
background of the ith data point. In general, the specific 
model used to correct for seeing and background systemat- 
ics is not important, provided that the model parameters 
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Fig. 3.— Top Panel: The cleaned PLANET data for OGLE-1998-BUL-14. As in Figure 1, the magnification A = (F — F B )/F S is plotted. 
Only PLANET data prior to HJD' = 1050, constituting ~ 95% of the complete PLANET data set are shown. The dashed line indicates 
the best-fit PSPL model, which has a time scale ts = 40 days, and a minimum impact parameter uq = 0.064, corresponding to a maximum 
magnification of j4 max ~ 16. The entire data set consists of 461 /-band and 139 V-band data points, the majority of which was taken between 
HJD' ~ 950 to HJD' ~ 1040. The median sampling interval during this time span is about 1 hour, or 10 — 3 £e, with no gaps greater than four 
days. The data are from the Yale-CTIO lm, South African Astronomical Observatory lm, the Perth 0.6m, the Canopus lm, and the CTIO 
0.9m. Bottom Panel: Residuals from the best-fit PSPL model. 



are not strongly correlated with rj and £. For most of the 
deviations we will be considering in this paper, i.e. those 
arising from nearly equal mass binary lenses, parallax, and 
finite source effects, model parameters are nearly uncorre- 
cted with r\ and £■ However, deviations caused by small 
mass ratio companions can occur over the course of several 
hours, and thus these deviations can be highly correlated 
with variations in seeing. Since these binary-lens models 
will not have the benefit of the additional fit parameters 
(77, £), the significance of any short-duration deviation will 
be overestimated. We will return to this point in § 5.2. 

We fit the model in the following way. We choose trial 
values for tE,uo, and to. This gives a prediction for the 
PSPL magnification as a function of time A(t). The two 
parameters Fg and Fb are then determined by performing 
a linear fit to the flux. This gives the PSPL flux mpsph(ti) ■ 
Finally, the parameters rj and £ are determined by per- 
forming a linear fit in magnitudes. The final x 2 f° r the 
trial model with parameters (£e, to, uq, i*s, Fb, rj, Q) is then 
evaluated in magnitudes and the values of these parame- 
ters that minimize \ 2 are then found using a downhill- 



simplex method ( Press ct al. 1992 ). Using the final val- 
ues of rj and Cj the data are corrected for the seeing and 
background systematics. This preliminary PSPL fit pro- 
duces x 2 ~ 2000 for 569 degrees-of- freedom (dof). Since 
no gross deviations from a PSPL light curve are appar- 
ent in the data, the high x 2 /dof is an indication that the 
DoPHOT errors are underestimating the true error, in this 
case by a factor of ~ 2. Some of this inflated \ 2 arises 
from a few highly-deviant outliers. Furthermore, different 
light curves have scatter that are underestimated by dif- 
ferent amounts. For these reasons, simply scaling all the 
errors by a factor that forces x 2 /dof to be unity would 
not be appropriate. We therefore adopt the following pro- 
cedure. Using the preliminary model, we first find the 
largest > 3a outlier and reject it. We recompute \ 2 f° r 
each light curve, and then rescale the errors in the individ- 
ual light curves by a factor that forces the x 2 /Ni for each 
light curve to be equal to unity, where Ni is the number 
of measurements in light curve i. All the errors are then 
scaled with an overall factor to force % 2 /dof for the en- 
tire data set to unity. This light curve is then refit to the 
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Fig. 4. — The histograms show the distributions of the residuals of the individual data points from the best-fit PSPL model, divided by 
their respective errors, before correction for seeing and background systematics (dashed histogram) and after correction (solid histogram). 
The solid curve is a Gaussian of unit variance. 



model in equation (Q), removing any residual seeing and 
background correlations. We iterate this process, succes- 
sively removing the largest 3cr outlier, rescaling the errors, 
and refitting the light curve, until no further 3<r outliers 
remain, and the fit has converged. The final scaling factors 
are given in Table 1. We find a total of five outliers that 
deviate from the final model by more than > 3cr. At the 
final step, we form two data sets. For the first data set, 
we reintroduce the > 3cr outliers, but with errors scaled 
by the same factor as their parent light curves; we will 
refer to this as the "cleaned" PLANET data set. The 
second data set does not contain the outliers; we refer to 
this as the "super-cleaned" PLANET data set. Finally, for 
some analyses, we combine the OGLE and super-cleaned 
PLANET data sets, scaling the OGLE errors by 2.01 so 
as to force x 2 /d.o.f. to unity. We will refer to this as the 
OGLE+PLANET super-cleaned data set. 

We now fit these corrected data sets to a PSPL 
model with parameters ts,, to, Uq, F$, Fb (Eq. [l]). For the 
PLANET datasets, there are 14 light curves, and thus the 
fit is a function of 31 parameters. Including OGLE data 
increases to the number of parameters by two. For the 
cleaned PLANET data set we find fit parameters and for- 
mal lcr errors determined from the linearized covariance 



matrix of t E = (39.6 ± 1.1) days, u = 0.0643 ± 0.0002, 
to = 956.016 ± 0.005. The fit parameters for the super- 
cleaned PLANET data set and the OGLE+PLANET 
super-cleaned data set are the same within the errors. 
The parameters and formal la errors for all three fits are 
summarized in Table 1, along with the blend fractions 
g = F B /F S for each light curve. The OGLE+PLANET 
dataset is shown in Figure 2, while the cleaned PLANET 
light curve and the best-fit PSPL model are shown in Fig- 
ure 3, along with the residuals from the model. 

In Figure 4 we show the distribution of residuals of the 
cleaned PLANET dataset from the best-fit PSPL model 
divided by their respective (rescaled) errors, along with a 
Gaussian of unit variance. Other than a few outliers, the 
two distributions are very similar, indicating that our er- 
rors are nearly Gaussian distributed. In order to illustrate 
the importance of including the seeing and background 
corrections, we also show the distribution of residuals di- 
vided by their respective rescaled errors before these cor- 
rections. The uncorrected data have a broader distribution 
with a median value that is systematically offset from zero; 
the entire distribution is highly non-Gaussian. 

In Figure 5 we show the residuals of the cleaned 
PLANET dataset from the PSPL model as a function of 
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Fig. 5. — The points with error bars are the residuals (in magnitudes) of the OGLE-1998-BUL-14 data points from the best-fit PSPL 
model as a function of the magnitude predicted by the model. We show the residuals separately for the /-band datasct (top panel) and the 
V-band dataset (bottom panel). The la scatter for the entire /-band data set, as denoted by the dashed line, is ~ 4%, whereas the scatter 
for all points with I < 16 is 1.5% (dotted line). The la scatter for the V-band data set is ~ 10% (entire data set), and ~ 7% (V < 18.7). 



the predicted magnitude of the model. The photometry in 
I is excellent; the la scatter for I < 16 is 1.5%; while the 
scatter for the entire / dataset is 4%. The scatter is ap- 
proximately 2.5 times larger in V; however this constitutes 
only 20% of the entire dataset. 

5. LIMITS ON COMPANIONS TO THE LENS 

The excellent coverage and high-quality photometry, 
combined with the fact that OGLE-1998-BUL-14 was a 
high- magnification event, make this an excellent candidate 
for the detection of planetary perturbations. Direct exam- 
ination of the light curve and residuals reveals no obvious 
planetary signatures, and, in fact, no obvious deviations 
from the PSPL model of any kind. However, since the 
deviations could be quite subtle, it is important that the 
light curve be searched systematically for any deviations. 
If no significant deviations are found, the good photome- 
try and excellent coverage can be used to place limits on 
the kinds of companions to the lens that could be present. 
To do this, we must calculate the detection efficiency of 
the OGLE-1998-BUL-14 light curve to companions. 



5.1. Detection Efficiency of OGLE-1998-BUL-U to 
Companions 

The efficiency e of a particular microlensing light curve 
to the detection of a binary system depends sensitively on 
two quantities: the the mass ratio of the system, q, and 
the angular separation d in units of 6b- The efficiency 
e(d, q) tends to zero when q — > 0, d — ► 0, or d — > oo (i.e. 
when the companion has a low mass compared to the pri- 
mary, or is very close or very far from the primary). We 
simultaneously search for planetary deviations and calcu- 
late the detection efficiency for OGLE-1998-BUL-14 using 
a method proposed by Gaudi & Sackett (2000). We briefly 
review the steps below. 

(1) The OGLE-1998-BUL-14 light curve is fit with a 

PSPL model by minimizing % 2 . The resulting % 2 
for this model is labelled Xpspl- 

(2) Holding d and q fixed, the binary lens model that 

best fits the observed light curve is found for 
each source trajectory a, leaving the 31 parameters 
(£e, uo, t 0} 14 x Fs, 14 x Fb) as free parameters. The 
difference Ax 2 (d, q, a) = x 2 {d, q, a) - Xp SPL is eval- 
uated. 
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Fig. 6. — The solid black lines are contours of constant detection efficiency, t{d,q), for all PLANET OGLE-1998-BUL-14 data, shown for 
projected separations d between the primary and companion in units of the Einstein ring radius, of < d < 4, and mass ratios between the 
primary and companion, q, of > log(q) > —5.0. The contours mark e = 5% (outer contour), 25%, 50%, 75%, and 95% (inner contour). 
The vertical dotted lines indicate the boundaries of the lensing zone, 0.6 < d < 1.6. The region of (d, q) parameter space shaded in black is 
excluded by the observations at the 95% confidence level. The panels are the results for different rejection criteria, Ax? hr h = 25 (top), 50 
(middle), 100 (bottom). A mass ratio of q = 10 — 2 corresponds to the mass ratio between a lOMj planet and a G dwarf. For this mass ratio, 
a companion with projected separation 0.4 < d < 2.4 is excluded at the 95% confidence level for all our rejection criterion. This corresponds 
to a range of instantaneous projected separations in physical units of 1.2 to 7.4 AU, assuming a G dwarf primary at 6.5 kpc. 



(3a) All parameter combinations (d, q, a) yielding 
A\ 2 (d, q, a) < A%g ag are nagged for further study as 
possible detections, where A%g ag is some reasonable 
detection criterion. 



(3b) The fraction of all binary-lens fits for the given (d, q) 
that satisfy the criterion A% 2 > Axth res h m com- 
puted, where Axthrcsh ^ s a rejection criterion. The 



detection efficiency e(d,q) of the data set for the as- 
sumed separation and mass ratio, 



1 f 27T 

e(d, q) = — da &[A X 2 (d, q, a) - Ax t 2 hrcsh ], 

(8) 

where Q[x] is a step function, is then computed. 
Note that a is uniformly distributed. 
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Fig. 7. — The detection efficiency of PLANET OGLE-1998-BUL-14 data to a companion averaged over the lensing zone, 0.6 < d < 1.6, 
as a function of the mass ratio q between the primary and the companion. We show the lensing zone detection efficiency for three different 
rejection thresholds, AXthrcsh = 25 (solid), 50 (dotted), and 100 (dashed). 



(4) Items (2) and (3) are repeated for a grid of (d, q) val- 

ues. This gives the detection efficiency e(d,q) for 
OGLE-1998-BUL-14 as a function of d and q, and 
also all binary-lens parameters (d, q, a) that give rise 
to significantly better fits to the light curve than the 
PSPL model. 

(5) The parameter combinations found in step (3a) are 

then used as initial guesses for our binary-lens x 2 
minimization routine, leaving all 34 parameters as 
free parameters, in order to find the minimum x 2 
and best-fit parameters for this local minimum. 

We search for binary-lens fits and calculate e(d,q) for 
< d < 4 and > log(q) > —5 at intervals of 0.1 in 
d and 0.5 in log(g). Since the search is performed on a 
grid of (d, q, a) and the grid points are unlikely to be sit- 
uated exactly on local minima, it important that |A%g ag | 
be sufficiently low so that all probable fits are found. We 
chose Ax| ag = —9 for our flagging criterion, and find only 

three combinations of d and q for which A% 2 < Axfl ag - 
Minimizing x 2 in the neighborhood of these trial solutions 
reveals that the best-fit parameters are quite close to the 
initial parameters, and % 2 decreases minimally. The ab- 



solute best-fit binary lens to OGLE-1998-BUL-14 has a 
X 2 = 720.6, or Ax 2 = —9.2. Since this is well below our 
threshold for detection of Ax 2 = —50, we do not claim 
a detection. Indeed, even a naive calculation based on 
Gaussian statistics and three additional parameters to de- 
scribe the companion yields a chance probability of 2.7% 
for Ax 2 < —9.2, which cannot be considered a detection 
at even the 3-sigma level. 

In fact, the probability of a random fluctuation of this 
magnitude is substantially higher than 2.7%. First, since 
we are including outlier points in the search, x 2 should be 
renormalized by the number of degrees of freedom, which 
would imply Ax 2 = —7.2 for which the chance probabil- 
ity with three additional parameters is 6.7%. However, 
the chance probability of a false planet detection is sub- 
stantially higher than this. Moreover, it cannot be com- 
puted from a x 2 table and would have to be determined 
by Monte Carlo simulation. To understand why, consider 
one of the planetary models we found that fits OGLE- 
1998-BUL-14 with Ax 2 = -7.2. The "success" of this 
model is basically driven by 15 points clustered within 
1 day whose mean lies ~ 0.6% above the PSPL model. 
The chance for such a deviation on this particular day is 
~ exp(Ax 2 /2)/(27r|Ax 2 |) 1/2 = 0.4%. Since there are a to- 
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tal of ~ 600 data points, there is a similar probability for 
such a fluctuation on each of 600/15 = 40 time intervals. 
Hence, the total probability is 1 - (1 - 0.004) 40 - 15%. 
In fact, the probability is somewhat higher still because 
we have taken account of only 15-point clusters and not 
larger or smaller clusters that could also mimic a planet. 

In any event, we have set our detection threshold sub- 
stanti ally higher than would be warranted solely to avoid 



is motivated by concern over unrecognized systematics 
which, experience has taught us, often give rise to spu- 
rious detections of formally high significance. We there- 
fore conclude that the light curve of OGLE-1998-BUL-14 
is consistent with a single lens within the uncertainties. 

5.2. Resulting Constraints on Companions 

The detection efficiency at a given (d, q) is the prob- 
ability that a companion of separation d and mass ratio 
q would have produced a deviation inconsistent (in the 
sense of A% 2 > Ax 2 hrcsh ) with the observed OGLE-1998- 
BUL-14 light curve. If no deviations are seen, compan- 
ions with (d, q) can be ruled out at a confidence level of 
e(d,q). The choice of the rejection threshold Ax 2 hresh 
can have a significant effect on the resulting detection 
efficiency, especially for low th resholds and mass ratios 
q< 10" 3 ( |Gaudi fc Sackett 2000| ). We choose to be conser- 
vative, and adopt Ax 2 hrosh = 50 as our fiducial threshold. 
For comparison, we also show the results for Ax 2 hrcsh 



thresh 

= 25 

and A Xt 2 hrosh 

The binary-lens detection efficiency of OGLE-1998- 
BUL-14 as a function of the mass ratio q and angular sepa- 
ration d is shown in Figure 6. The darkest shading denotes 
those parameter combinations (d, q) for which e > 95%, 
and thus are excluded at the 95% confidence level from 
lying above the rejection criterion. Table 2 shows the 
range of angular separations d that are excluded by our 
OGLE-1998-BUL-14 data set for several mass ratios and 
the three different rejection thresholds, Ax 2 hrcsh = 25,50, 
and 100. For q> 10" 15 , e = 100% out to the largest sepa- 
ration for which we calculate e, namely d = 4.0. For these 
mass ratios d = 4 is thus a lower limit to the excluded 
range; the upper end of the excluded range is likely to be 
considerably larger. We find that any companion to the 
primary lens with mass ratio q~ 10~ 2 and angular separa- 
tion 0.4 < d < 2.4 is excluded by the data, that is, such 
a companion would produce deviations at the A\ 2 > 100 
level that are not observed. 

Note that these limits apply to individual companions 
only, not to systems of companions. Implicit in our cal- 
culation of e(d, q) is the assumption that multiple planets 
affect the magnification structure of the lens, and there- 
fore the deviation from the single-lens magnification, in 
an independent way. This assumption is likely to break 
down in regions near the central caustic when more than 
one p lanet is in the lensing zone (Gaudi, Naber, & Sackett 
1998). In this case, the efficiencies calculated using the 
method of Gaudi & Sackett (2000) will be in error by an 
amount that will depend on the relative mass ratios, orien- 
tations, and projected separations of the two companions. 
Since OGLE-1998-BUL-14 is a high-magnification event, 
and most of the constraints come from portions of the 
light curve near the peak of the event (i.e. where the cen- 
tral caustic is probed), our results are only strictly valid 



for single planets. A full investigation of the effect of mul- 
tiple planets of e is beyond the scope of this paper. We 
expect, however, that the planet detection efficiencies for 
multiple Jovian planetary systems may actually be higher 
than for single planets near the peak of OGLE-1998-BUL- 
14, since the region of anomalous magnification near the 
central caustic generally occupies a larger fraction of the 



chanc 3 statistical fluctuations. This conservative approach & Sackett 1998) 



Einstein ring wh en two planets are present (Gaudi, Naber 



It is interesting to compare the limits on companions 
to the primary lens of OGLE-1998-BUL-14 to those for 
the primary lens of MACH O-98-BLG-35 der ived by the 
MPS/MOA collaborations ( |Rhic et al. 1999b| ). MACHO- 
98-BLG-35 was a higher magnification event (A max ~ 75) 
than OGLE-1998-BUL-14 (,4 max ~ 16). Since higher mag- 
nifica tion events have higher i n trinsic detection efficien - 
cies ( priest fc Safizadch 199^ ; paudi fc Sackett 2000| ), 
one would expect, for similar sampling and photometric 
precision, the limits on companions to be more stringent 
for MACHO-98-BLG-35. However, although the photo- 
metric precision obtained by MPS/MOA on MACHO-98- 
BLG-35 is similar to that for OGLE-1998-BUL-14 (~ few 
percent), the sampling of MACHO-98-BLG-35 (in terms 
of fraction of the time scale £e) is poorer, due primarily 
to the fact that MACHO-98-BLG-35 was a shorter time 
scale event. Although Rhie et al. (1999b) used a slightly 
different method to calculate e(d, q) than that suggested 
by Gaudi & Sackett (2000), and used a different rejection 
criterion (Ax 2 h ros h = 40), we can make a rough compar- 
ison between their resulting detection efficiencies shown 
in their Fig. 7 with ours for Ax 2 hrosh = 50 shown in the 
middle panel of Fig. 6. We see that the detection efficien- 
cies to companions are everywhere higher for MACHO- 
98-BLG-35 than for OGLE-1998-BUL-14. This indicates 
that, when the peak of the event can be measured, maxi- 
mum magnification is a more important factor than sam- 
pling in determining the constraining power (and hence the 
ability to detect companions) in an observed microlensing 
event. 

In Figure 7 we show the detection efficiency averaged 
over the lensing zone (where the detection efficiency is the 
highest), 0.6 < d < 1.6, 



£Lz(g) 



1.6 



e(c?, q) dd 



(9) 



for mass ratios < q < 10~ 5 . The lensing zone detection 
efficiencies for several representative mass ratios are tabu- 
lated in Table 3. For a planetary model in which com- 
panions have angular separations distributed uniformly 
throughout the lensing zone, clz represents the probability 
that a companion of mass ratio q would have been detected 
with the OGLE-1998-BUL-14 data set. These probabili- 
ties are quite high: for example, the detection efficiency 
for a companion of mass ratio q~ 10 -3 in the lensing zone 
of the OGLE-1998-BUL-14 primary is 2 80%. 

For this analysis, we have assumed that the source can 
be treated as point-like. Finite source effects will have a 
substantial effect on e(d, q) if the angular size of source 
= 9*/9-e is comparable to the Einstein ring radius of 
the companion, 9 p = B^q 1 ^ 2 ( Gaudi fc Sackett 2000 ). Any 



planetary deviations will be broadened but reduced in am- 
plitude for p^q 1 ! 2 . For OGLE-1998-BUL-14, no finite 
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source size effects were detected, thus we can piace an up- 
per limit on p*. As we show in the § 6.2, the 3c limit 
is < 0.062 . Thus the detection efficiencies calculated 
above are strictly valid only for q~ 10 -2 5 . However, for 
typical lens parameters, is likely to be considerably 
smaller, ~ 0.01. Statistically, we thus expect the re- 
sults to be valid for q~ 10 . For mass ratios less than 
this, the detection efficiencies calculating using a point 
sourc e may be overestimated by tens of percent (Gaud] 
& Sa dkett 20001 ). 

Poorly-constrained blend fractions can also induce sub- 
stantial uncertainties in the derived detection efficiencies, 
due to the correlation between bl ending and impact pa- 
rameter wo ( paudi fc Sackctt 200C ) . Fortunately, blending 
is easier to constrain in high-magnification events. Indeed, 
due to the high magnification of OGLE-1998-BUL-14 and 
the dense and precise photometry, the blend fraction of 
the event, and therefore uo, are quite well constrained. 
The fractional error in uq is < 1% and should contribute 
negligibly to the uncertainty in e(d, q). 

When fitting for the binary-lens model, we did not in- 
clude seeing and background correlation terms (§ 4). Al- 
though the proper approach would be to include these 
terms in all binary-lens fits (and indeed, all fits in gen- 
eral), the computational cost is prohibitive. Since these 
terms were included in the PSPL fit, this implies that 
A\ 2 between the binary-lens fit and single-lens fit may 
be overestimated if the deviations from the PSPL fit due 
to the binary lens are highly correlated with either the 
seeing or background. This is generally not a problem for 
background correlations because the systematic deviations 
arising from the long time scale changes in the background 
are unlikely to be confused with deviations caused by com- 
panions. However seeing correlations are more insidious, 
because the short-time scale deviations caused by low mass 
ratio companions could be confused with systematic devi- 
ations arising from the nightly seeing changes in a poorly 
sampled light curve. Consider a deviation from the PSPL 
light curve that is perfectly correlated with the seeing. It 
can be shown that the fractional error in \ 2 one makes by 
not including the seeing correlation in the binary lens fit 
is A\ 2 /x 2 ~ Ns/Ntot, where iV^ is the number of deviant 
data points and A to t is the total number of data points 
in the light curve. For uniformly sampled data, this is 
simply the ratio of the time scale of the perturbation t$ 
to the total duration of observations. For planetary mi- 
crolensing, t$ ~ q 1 / 2 ^, and thus Ax 2 /x 2 ~ q 1 ^ 2 - Thus 
for small mass ratios, q < 10~ 3 , the fractional error in x 2 
is small ~ 3%. For q > 10~ 3 , the companions produce co- 
herent deviations that last many days, and thus cannot be 
correlated with the seeing. We therefore conclude that ne- 
glecting the seeing and background correlations does not 
result in seriously overestimated detection efficiencies. 

6. LIMITS ON THE MASS OF THE LENS 

In the previous section, we placed limits on possible 
companions to the primary lens responsible for the mi- 
crolensing event OGLE-1998-BUL-14 by calculating the 
detection efficiency as a function of the mass ratio q and 
angular separation d between the primary and secondary 
in units of the angular Einstein ring of the system. Al- 
though limits on the mass ratio are interesting in their 
own right, what is of interest ultimately are limits on the 



mass of the companion, M p , and the physical orbital sep- 
aration of the companion, a. In order to translate limits 
on q and d into limits on M p and a, one must know the 
mass of the primary lens, M, and its distance Dl, and the 
orbital phase and inclination of the system. 

Unfortunately it is not generally possible to measure the 
mass and distance to the lens. The only parameter one can 
measure from a generic single lens event that contains in- 
formation about the lens is the time scale £e, a degenerate 
combination of the mass M, distance D^, and transverse 
velocity v of the lens (c.f. Eq. ||). However, detection 
of various higher order effects, such as parallax or finite 
source effects, enables one to extract additional informa- 
tion and partially break this degeneracy. In the following 
subsections, we search for signatures of these effects in the 
light curve of OGLE-1998-BUL-14. Other than a marginal 
detection of a parallax asymmetry, we do not detect either 
of these effects, and use their absence to exclude regions 
of the (M, £>l) plane. In addition, we place a limit on the 
amount of light emitted from the lens itself, and translate 
this into a limit on the lens mass. As we show, despite the 
excellent photometry and coverage of OGLE-1998-BUL- 
14, the limits on these higher order effects are not very 
stringent, and do not translate into strong constraints on 
the mass of the lens. 

The analysis uses the super-cleaned PLANET data set 
since, as opposed to planetary perturbations, the effects we 
are searching for here induce subtle deviations with time 
scales of many days, and therefore cannot be described by 
a isolated large outliers. Using the super-cleaned data set 
(which does not include these outliers) will provide more 
robust and reliable limits. For the parallax analysis only 
we use the OGLE+PLANET cleaned data set in order to 
have a more secure description of the total baseline flux. 

6.1. Parallax Limits 

The motion of the Earth around the Sun induces de- 
partures from rectilinear motion in the path of the lens 
relative to the Earth-source line-of-sight and thus alters 
the formula for the angular separation u(t) as a function 
of time. This phenomenon is commonly referred to as par- 
allax, and in general gives rise t o light curves that deviate 
from the standard PSPL model ( Gould 1992| ; Alcock et al. 
1995 ). The magnitude of the deviation depends on the 
time of year when the event peaks, the angle tfi of the lens 
trajectory with respect to the north ecliptic pole, and the 
transverse velocity v = v(Z?s/-Dls) of the lens projected 
onto the observer plane. 

If the time scale of the event is short compared to the 
period of the Earth's orbit, the Earth's acceleration can be 
approximated as constant over the course of the e vent. In 
this case, the angular separation c an be written as ( |Gould 
Miralda-Escude, fc Bahcall 1994|), 



i(t) = [e+ul] 1/2 , ^t+^ e t 2 , (10) 



where r = (t — to)/t-E. The asymmetry parameter n is 
given by, 

(11) 



k = Q ffi — — sin A sin < 



where Q§ = 2ir yr , v§ ~ 30 kms 1 is the linear speed 
of the Earth around the Sun, and A is the angle between 
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Fig. 8. — Excluded regions for the lens mass and distance based on the lack of second-order signatures in the light curve of OGLE-1998- 
BUL-14. The solid lines show 3<r upper or lower limits to the mass of the lens in solar masses M/Mq as a function of the distance of the 
lens in kpc for a source distance of 8 kpc. The solid hatching shows the excluded region from the lack of parallax effects. The short-dashed 
hatching shows the region excluded by the lack of finite source effects. Combining the finite source and parallax limits, we obtain a lower 
limit to the mass, which is shown as the triangle. The long-dashed hatching shows the region excluded if the lens is a main-sequence star; the 
upper limit to the lens mass is shown as a cross. The dotted horizontal line is the hydrogen-burning limit. Our adopted lens mass M = Mq 
and distance Z?l = 8 kpc are shown as the square. 



the source and Sun at the time of maximum magnification. 
In the case of OGLE-1998-BUL-14, sin A - 0.5. Implicit 
in equation ( pT| ) is the approximation that the source is in 
the Galactic plane. Thus, for short-time scale events one 
can measure only the degenerate combination v/ sin cf>. 

We fit our OGLE-1998-BUL-14 light curve data with 
the parallax asymmetry model above and find an improve- 
ment of Ax 2 = 4.4 for 1 extra dof over the standard PSPL 
model, corresponding to a ~ 2a detection of an asymme- 
try. We find a best-fit value of re = 6.1 ± 3.6 x 1CT 4 (ler), 
corresponding to 



(420 ± 250) km s" 1 (la). 



(12) 



This detection is not particularly useful for our purposes, 
however, not only because it is merely a 2a detection, but 
also because the direction of lens motion is unknown. 

In order to obtain a constraint on v that is indepen- 
dent of <fi, we must use the form for u(t) that inclu des 
the full two-dimensional (v,cj)) parallax information (Al- 



cock et al. 1995). Given that the detection of asymmetry 



was marginal, it is not surprising that we are unable to 
obtain independent constraints on v and <fi. Calculating 
X 2 as a function of v, and letting all other parameters (in- 
cluding (f>) vary, we recover the detection of the parallax 
asymmetry for v~ 30 kms -1 . For v~ 30 kms -1 , however, 
v becomes comparable to the Earth's velocity. For these 
small velocities, the component of the Earth's velocity per- 
pendicular to the motion of the lens becomes significant, 
and the light curve begins to deviate appreciably from the 
PSPL model, in a manner that is inconsistent with the ob- 
servations for all values of <f>. This produces a lower limit 
to the projected velocity that is independent of <f>, 



> 28.5 kms" 



(3a). 



(13) 



We can combine this limit with the time scale of the event 
to give a lower limit to the lens mass as a function of the 
relative lens distance, 



M(x) > 6.7 x 10 -3 M Q 



1 



(Parallax, 3a) (14) 
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Fig. 9. — Calibrated, dereddened color-magnitude diagram for the field centered on OGLE-1998-BUL-14 from Yale-CTIO lm data. The 
mean reddening for the field is Aj and E(V — 7)=1.56, determined from the mean / magnitude and V — I color of the clump. The arrow 
shows the direction and magnitude of the dereddening vector. The open circle shows the dereddened position of the unmagnified source, 
which is composed of the microlensed source (filled circle) and an unresolved blended star (cross). The dotted line shows the 3<r upper limit 
to the dereddened I magnitude of the lens, assuming the extinction to the lens is the same as the OGLE-1998-BUL-14 field. 



where x = Z?l/ D§. This limit is shown in Figure 8. 

For bulge self-lensing, both lens and source belong to a 
population with approximately isotropic velocity distribu- 
tions, producing no preferred direction for the transverse 
velocity, v, and thus no preferred value for (j>. For disk 
lenses, however, there is a preferred direction for <f> due to 
Galactic rotation, in which case we would obtain a stronger 
limit 



. Uufui luiialely, we do uoL know apiiuii whether 
the le hs ia in the bulge or disk. We will therefore adopt 

the conservative assumption that there is no preferred di- 
rection for v, and use the limit given in equation (|l4|). 



6.2. Finite Source Limit 

A point-lens transiting the face of a source will resolve 
it, creating a distortion in the magnification that deviates 
from the form given in equation (|2]). A detection of this 
distortion gives a measurement of the angular size of the 
source, f?*, in units of the angular Einstein ring radius, 
p„ = 6J9 E ([Gould 1994| ; [Ncmiroff k Wickramasinghc 
1994| |Witt fc Mao 1994|). The requirement for such finite 



source effects to be detectable is that the impact parame- 
ter of the event must be comparable to or smaller than the 
dimcnsionlcss source size, uq~P*- No such finite-source 
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deviations are apparent in the light curve of OGLE-1998- 
BUL-14. This implies an upper limit to the dimensionless 
source size of p* ~ 0.06. For a more exact limit, we calcu- 
late x 2 as a function of p*, leaving all other parameters 
free to vary, and find that 



p* < 0.062 (3a). 



(15) 



In calculating this limit, we have assumed a uniform sur- 
face brightness profile for the source. Using a more realis- 
tic, but model-dependent, limb-darkened profile weakens 
this limit slightly. 

In order to convert this upper limit on p* into a lower 
limit on $e, we must know the angular size of the source, 
6**, which can be estimated from its (V — I) color and 
/ magnitude. Figure 9 shows the calibrated dereddened 
color magnitude diagram (CMD) for OGLE-1998-BUL-14 
field from Yale-CTIO lm data. The stars in the OGLE- 
1998-BUL-14 field were calibrated relative to observations 
of 9 standard stars in Landolt (1992), measured several 
times during the night of 14 August 1998 at the Yale- 
CTIO lm telescope in the same V and / filters used for 
the observations of the microlensing event. Extinction and 
color correction terms were derived in the normal manner. 
Comparison with similar measurements on other nights 
suggests the accuracy of the calibration is ~ 0.02 mag. 
We fit the distribution of I magnitudes and (V — I) colors 
of the observed clump to the model of Stanek (1995). We 
then determined the E(V — I) and Aj by comparing our 
fitted clump magnitude and color to the dereddened / c i,o 
and (V — I) c \,o °f the clump for the bulge as determined by 
Paczyhski &Stanek (1998). We find E(V - I) = 1.56 and 
Aj = 2.16. Using the PSPL fit and Yale-CTIO standards, 
we determine the calibrated, dereddened color and mag- 
nitude of the microlensed source to be Iq = 14.73 ± 0.03 
and (V — I)q = 1.23 ± 0.08, where the errors reflect the 
calibration and model uncertainties added in quadrature. 
From its position on the CMD, Figure 9, we conclude that 
OGLE-1998-BUL-14 is likely to be a clump giant or a RGB 
star. 

To obtain an estimate of the source radius, 0*, we use 
the empirical color-surface brightness relationship for gi- 
ant stars derived by van Belle (1999), which we rewrite 
as 

0, = 6.0 ^as 10 -o.2 ( v-i6)+o.5[(v-/)-i.2] : (16) 

where we have assumed the relationship (V—K) = 2.2(V— 
I), derived from Bertelli et al. (1994) isochrones. For our 
parameters, we find = 6.3/ias. An uncertainty in the 
extinction S Ay leads to an uncertainty in the angular size 
of 59*/ 8* ~ —0.16 5 Ay, where the coefficient dep ends on 
the temperature of the source ( Albrow et al. 1999 ), which 
we have assumed to be T — 5000 K. We estimate the 
uncertainty in the extinction to be S Ay ~ 0.15, based on 
the dispersion of the clump, resulting in an uncertainty 
of ~ 2% in 9*. The 3% uncertainty in the /-magnitude 
of the source leads to an additional uncertainty of ~ 1% 
in 9*; the 8% uncertainty in the color leads to an uncer- 
tainty of ~ 6% in 9*. Adding these in quadrature, we find 
9* = (6.3 ± 0.4) ^as. 

Adopting this value for the angular size of the source, 
we can now translate the limit on p* (Eq. 15) directly to 
a limit on the angular Einstein ring radius, 



Since #e depends only on the lens mass and distance and 
on the distance to the source, this limit can be written 
in terms of a limit of the lens mass as a function of the 
relative lens distance, 



M(x) > 1.0 x 1O~ 2 M 



1 



(Finite Source, 3<r) (18) 



We combine theparallax limit (Eq. |l4|) and the finite 
source limit (Eq. [18|) to obtain a lower limit to the lens 
mass M > 8.3 x l(U 5 Af Q , which occurs at x = 3.55. This 
lower limit is indicated in Figure 8. 

6.3. Luminous Lens Limit 

If the lens of OGLE-1998-BUL-14 is a main-sequence 
star, it emits light. Although the angular separation be- 
tween the lens and source is much too small for the lens to 
be resolved, additional light from the lens could, in princi- 
ple, alter the shape of the light curve (c.f. Eq. [l]). For most 
microlensing events, degeneracies between the fit parame- 
ters make it difficult to constrain accurately the amount of 
blended light, and thus any light that may be arising from 
the lens itself. Fortunately, the good photometry, complete 
coverage, and small uq of OGLE-1998-BUL-14 enable the 
blend fraction to be constrained quite tightly, allowing us 
to place an interesting upper limit on the blended light 
emitted by the lens. 

To do this, we rewrite the flux of any unresolved light 
as Fb = Fu + Fl, where Fl is the flux of the lens, and Fj 
is the flux of any additional unresolved source. Equation 
ll is then 







F(t) = F s A(t) + Fu + F L . 



(19) 



We assume that the /-band flux of the lens is the same for 
all the /-band light curves, and similarly that the U-band 
flux is the same for all the U-band light curves. This is a 
reasonable assumption, since all the observatories use the 
same / and V filters, and because the lens and source will 
be unresolved for all observatories. However, due to dif- 
ferences in pixel size and image quality, we must allow ad- 
ditional unresolved light, Fu not associated with the lens 
or source to vary from site to site. We then compute x 2 
as a function of Fl, allowing the other parameters to vary, 
and imposing the constraint that Fb always be positive. 
For Fl < fsi Ax 2 between the fit with and without a lu- 
minous lens is small. However, as Fl increases, Ax 2 rises 
because Fj must be negative to match the observations, 
which is unphysical. In this way we find the minimum 
/-band magnitude of the lens that is consistent with our 
observations to be 



h > 



21.69- (V -I) 
19.94 



if (V - /) < 1.75 
if (V — I)> 1.75 



(20) 



6> E > 100 ^as (3a). 



(17) 



where the limit is 3a. Note that this limit is not dered- 
dened since the lens may have less extinction than the 
mean extinction toward the OGLE-1998-BUL-14 field. For 
reference, the dereddened limit is shown in Figure 9, as- 
suming the reddening to the lens is the same as for the 
OGLE-1998-BUL-14 field. 

In order to convert this limit on lens light into a limit 
on the lens mass as a function of its distance, we must 
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adopt a relationship between mass and / and V magni- 
tude. For this purpose, we use the solar metallicity Bertelli 
et al. (1994) theoretical isochrone s for M > 0.6 Mr, the 
solar metallicity Yale isochrones (Demarquc ct al. 1996) 
for 0.35 M Q < M < 0.6 M Q , and extrapolate the Yale 
isochrones using a 2nd order polynomial for < M < 
0.35 M Q . We assume that the dust is distributed uni- 
formly between the observer and 1.5 kpc, with a total 
reddening at 1.5 kpc equal to the mean reddening of the 
OGLE-1998-BUL-14 field, E(V-I) = 1.56 and Aj = 2.16. 
For all distances between and 8 kpc, we find the largest 
mass that is consistent with the apparent magnitude limit 
(Eq. |2(]). This mass-distance limit is shown in Figure 8. 
It can be adequately represented by, 

log[M(a:)] < 0.8+0.07a;-0.75:zr - 2 (Luminous Lens, 3er) 

.( 21 ) 

Varying the age, metallicity, or dust distribution within 
reasonable limits changes the mass limit at any distance 
by Z 0.2 dex. 

If the lens is a main-sequence star, the largest mass con- 
sistent with the lack of appreciable light in the light curve 
of OGLE-1998-BUL-14 is M ~ 1.3 M which occurs when 
£>l ~ 8 kpc (x = 1). Thus the lens must be a G dwarf or 
later. Of course, the lens may not be a main sequence star. 
Gould (1999) estimates that ~ 80% of events detected 
toward the bulge are due to main-sequence lenses. The 
remaining 20% are due to stellar remnants, i.e. as white 
dwarfs, neutron stars, or black holes. If the lens is such a 
remnant, it would not obey the mass- luminosity relation- 
ship used to find the limit in equation (^l|). However, white 
dwarfs and neutron stars have masses of Mwr> ~ 0.6 Mq 
and Mns ~ 1.35 M , and so automatically satisfy this 
limit for nearly all likely distances to the lens. Thus the 
vast majority (99%) of lenses will satisfy the limit in equa- 
tion @. 

6.4. Combined Limits 

The limits on the mass and distance to the lens set in 
the previous subsections are, for the most part, model in- 
dependent. Unfortunately, they are also not very strin- 
gent. Even with the assumption that the lens is a main- 
sequence star, the allowed regions in the (M, Di) plane 
are quite large, spanning two orders of magnitude in mass, 
0.01 Mq ~ M ~ 1.3 Mq, and nearly the entire range in dis- 
tance, 0.5 kpc~£>L~8 kpc. Our analysis indicates that, 
even with excellent coverage and good photometry, it will 
be quite difficult to routinely obtain stringent limits on 
the mass and distance to the lens for most events based 
on photometry alone. 

It has been shown by Dominik (1998) how probability 
densities for physical quantities of the lens system can be 
derived under the assumption of statistical distributions of 
the mass spectrum, the mass density, and the transverse 
velocity. Rather than doing this, we will simply note that 
if the lens is in the bulge (6 kpc~Z?L~8 kpc), and has 
a typical transverse velocity for bulge self-lensing events 
(v ~ 100 kms -1 ), then the measured £e implies that it 
is likely to have a mass near the upper end of the allowed 
range. However, we cannot rule out that the lens is moving 
slowly, and therefore that the mass is quite small. 

7. DISCUSSION 



7.1. The Detection Efficiency as a Function of Mass and 

Separation 

In order to convert the limits on companions in the q — d 
plane to limits on companions in the M p — a plane, we need 
estimates of the mass and distance to the lens. However, 
as we demonstrated in § 6, it is quite difficult to obtain 
stringent limits on these quantities from photometric data 
alone. For illustrative purposes, therefore, we will simply 
assume that the lens is a G dwarf and adopt M = M Q , and 
a lens distance of 6.5 kpc, so that te = 3.1 AU. We stress, 
however, that this choice is somewhat arbitrary, and that 
the lens mass may be smaller by two orders of magnitude. 

Since microlensing is only sensitive to the instantaneous 
angular separation, d, we must first convolve the detection 
efficiency e(d) with the distribution of d for a given semi- 
major axis a. To do this we integrate over all random incli- 
nations and orbital phases, assuming circular orbits. This 
distribution is given explicitly in Gould & Loeb (1992). 
Convolving the resulting distribution with e(d) gives the 
detection efficiency as a function of mass ratio and physi- 
cal (three dimensional) separation in units of j-e. We then 
use the values of M and te above to convert to e(q, a/r^) 
to the desired detection efficiency e(M p , a) as a function 
of the mass and true orbital separation of the companion 
in AU. 

This detection efficiency e(M p ,a) as a function of phys- 
ical parameters for our assumed primary lens (M = Mq 
and Z?l = 6.5 kpc) is shown in Figure 10, for our fiducial 
rejection threshold of Axt hrcsh = 50. Stellar companions 
to the primary lens of OGLE-1998-BUL-14 with separa- 
tions between ~ 2 AU and 11 AU, (the largest separa- 
tion for which we calculate e) are excluded. Although we 
cannot exclude a Jupiter-mass companion at any separa- 
tion, we have a ~ 80% chance of detecting such a com- 
panion at 3 AU. The detection efficiency for OGLE-1998- 
BUL-14 is > 25% at a = 3 AU for all companion masses 
M p > 0.03 Mj. We find that we had a ~ 60% chance of 
detecting a companion with the mass and separation of 
Jupiter (M p = Mj and a — 5.2 AU), and a ~ 5% chance 
of detecting a companion with the mass and separation of 
Saturn (M p = 0.3 Mj and a = 9.5 AU) in the light curve 
of OGLE-1998-BUL-14. 

Thus, although Jupiter analogs cannot be ruled out in 
OGLE-1998-BUL-14, the detection efficiencies are high 
enough that future non-detections in several events with 
similar quality will be sufficient to place meaningful con- 
straints on their abundance. 

7.2. Comparison with Other Methods 

How do the OGLE-1998-BUL-14 efficiencies compare to 
planet detection via other methods? In Figure 10 we show 
the radial velocity detection limit on M p sin i for a solar 
mass primary as a function of the semi-major axis for a ve- 
locity amplitude of K = 20 m s" 1 , which is the limit found 
for the majority of the stars in th e Lick Planet Search 
( Cummings, Marcy, fc Butler 199^ ). Although we show 
this limit for the full range of a, in reality the detection 
sensitivity extends only to a ~ 5 AU due to the finite du- 
ration of radial-velocity planet searches and the fact that 
one needs to observe a significant fraction of an orbital pe- 
riod. The limit rises dramatically for a~ 5 AU because the 
period of the companion becomes larger than the dura- 
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Separation (AU) 



Fig. 10. — Example detection efficiencies for the PLANET data set of OGLE-1998-BUL-14 as a function of the mass and orbital separation 
of the companion assuming a primary lens mass of Mq and Einstein ring radius of T'e = 3.1 AU. The contours are the same as in Figure 6. 
In order to convert from mass ratio and projected separation to mass and physical separation, we have averaged over orbital phase and 
inclination (assuming circular orbits). Jupiter and Saturn are marked with stars, as are the extrasolar planets discovered with radial velocity 
techniques. The horizontal line marks the hydrogen-burning limit. The dotted line shows the radial-velocity detection limit for an accuracy 
of 20 ms~ 1 and a primary mass of Mq. The dashed line is the astrometric detection limit for an accuracy of 1 mas and a primary of mass 
Mq at 10 pc. 



tion of the observations. In addition, we plot in Figure 10 
the M p sin i and a for planetary candidates detected in 
the Lick survey. Radial velocity searches clearly probe a 
different region of parameter space than microlensing, in 
particular, smaller separations. Note, however, that our 
OGLE-1998-BUL-14 data set gives us a > 75% chance of 
detecting analogs to two of these extrasolar planets: the 
third companion to Upsilon And and the companion to 14 
Her. Although the efficiency is low, we do have sensitivity 



to planets with masses as small as ~ O.OlAfj, considerably 
smaller than can be detected via radial velocity methods. 

For comparison, we also show in Figure 10 the astromet- 
ric detection limit on M p for a Mq primary at 10 pc, for an 
astrometric accuracy of ca = 1 mas. For an astrometric 
campaign of 11 years, this limit extends to ~ 5 AU. Such 
an astrometric campaign (a a = 1 mas, P — 11 years), 
would be sensitive to companions similar to those excluded 
in our analysis of OGLE-1998-BUL-14. The proposed 
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Space Interferometry Mission (SIM) promises ~ 4 /j,as as- 
trometric accuracy, which would permit the detection of 
considerably smaller mass companions. 

8. CONCLUSION 

We have presented the PLANET photometric data set, 
consisting of 461 /-band and 139 U-band measurements, 
for the microlensing event OGLE-1998-BUL-14. The me- 
dian sampling interval of one hour, RMS scatter of 1.5% 
over the peak, and high magnification of OGLE-1998- 
BUL-14 make this data set especially sensitive to the pres- 
ence of lensing companions. Within our photometric un- 
certainties, the data set is consistent with a single lens. 

Our analysis indicates that no companions with mass 
ratios q > 0.01 and instantaneous projected separations 
0.4 < d < 2.4 are present. Assuming a solar-mass primary, 
this mass ratio range includes known stellar binaries and 
super Jupiters (M p = 10 Mj). Less massive companions 
and those at larger or smaller separations are excluded 
with less confidence. Massive companions with q~ 10~ 15 
can be excluded for projected separations at least as large 
as 4 te- 

The absence of strong parallax, proper motion, or lens 
light detections allows us to constrain the mass of the lens 
to 0.01 Mq -MS 1.3 M©. Assuming a solar-type lens with 
M = Mq at a distance Z?l = 6.5 kpc, the Einstein ring 
radius of the primary corresponds to te = 3.1 AU. Us- 
ing this value, the PLANET light curve has efficiencies 
of 60% and 5% for Jupiter and Saturn analogs, respec- 
tively, and a greater than > 75% efficiency for compan- 
ions like those in the Upsilon And and 14 Her systems. 
Planets with M p > 10 Mj and true orbital separations 
1.2 AU < a < 7.4 AU are excluded, assuming these fidu- 
cial primary lens parameters. 

In performing our analysis, we have considered the sys- 
tematic effects of (1) correlations between our photome- 
try and image quality and sky background; (2) underes- 
timated error bars; (3) finite source size; and (4) poorly- 
constrained blending. We find that the DoPHOT-reported 
uncertainties underestimate the true scatter, and that the 
residuals from the best-fit model are significantly corre- 
lated with image quality and background. Applying a sim- 
ple linear correction term, rescaling the uncertainties, and 
eliminating five non-sequential outliers, however, results 
in a Gaussian error distribution yielding x 2 /dof = 1 for a 
point-source point-lens fit. This allowed us to proceed with 
a % 2 -analysis to find the portion of projected separation- 



mass ratio (d-q) binary parameter space excluded by our 
data. We estimate that non-zero source size is unlikely to 
significantly affect our results for qZ 10~ 4 and almost cer- 
tainly not for q~ 10~ 25 . Blending (and thus the true im- 
pact parameter) is well constrained for OGLE-1998-BUL- 
14 by our data set, and thus has a negligible effect on our 
conclusions. 

Our data set for this microlensing event is sensitive to 
planets occupying a different range of parameter space 
than current planet searches by other techniques. In par- 
ticular, super- Jupiters (M p ~ 10 Mj) can be ruled out as 
companions to the stellar-mass lens of OGLE-1998-BUL- 
14 at distances of several AU, larger than those probed 
by more sensitive - but relatively recently commenced - 
radial velocity and astrometric searches. 

It is not possible to derive general inferences about the 
abundance and characteristics of binary or planetary sys- 
tems from observations of any single system. Nevertheless, 
our results for OGLE-1998-BUL-14 clearly demonstrate 
the ability of microlensing to contribute to our knowledge 
of Jovian planets several AU from their parent stars, and 
- if data of high enough quality can be collected for a large 
enough number of events - to the search for and study of 
planets of much smaller mass as well. The analysis pre- 
sented here for OGLE-1998-BUL-14 represents the first 
step in the larger task of performing a combined analy- 
sis of the growing PLANET data base of frequently and 
precisely monitored microlensing light curves. When com- 
pleted, statistical inferences can be drawn about the fre- 
quency and distribution of stellar and Jovian companions 
to stellar lenses in the Galaxy. 
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real-time alerts and for making their data publicly avail- 
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cluding Andreas Berlind, Alberto Conti, Paul Martini, 
and Andrew Stephens. This work was supported in part 
by NASA under Award No. NAG5-7589, the NSF un- 
der grants AST 97-27520 and AST 95-30619, Nederlands 
Wetenschapelijk Onderzoek through award GBE 614-21- 
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Filter 


Site 


Number of Points 


o-/o- D oP a 


/ 


CTIO 0.9m 


48 


1.55 




Perth 


50 


1.00 




Canopus A 


60 


1.78 




Canopus B 


47 


2.25 




Yale 


56 


1.83 




SAAO A 


113 


1.43 




SAAO B 


20 


1.43 




SAAO C 


67 


1.42 




Total PLANET I 


461 


— 


V 


Perth 


4 


0.08 




Canopus A 


20 


1.77 




i ale 


00 


l.OO 




ca An A 


oo 


1 .Oo 




SAAO B 


6 


0.94 




SAAO C 


20 


0.88 




Total PLANET V 


139 






Total PLANET I+V 


600 




I 


OGLE 


159 


2.01 




Total PLANET+OGLE 


754 





Table 1 

Number of Data Points and Error Scaling Factors 



"The scaling factor for the DoPHOT reported errors. 







PLANET CL a 


PLANET SC b 


OGLE+PL SC C 


to 




956.016±0.005 


956.011±0.005 


956.011±0.005 


*E 




39.6±1.1 


40.0±1.2 


40.0±0.58 


u 




0.0643±0.0002 


0.0639±0.0002 


0.0639±0.0002 


9i 


CTIO 0.9m 


0.16±0.05 


0.17±0.05 


0.17±0.04 




Perth 


0.00±0.03 


0.00±0.03 


0.00±0.02 




Canopus A 


0.03±0.04 


0.04±0.04 


0.04±0.04 




Canopus B 


0.14±0.05 


0.14±0.05 


0.14±0.04 




Yale 


0.03±0.03 


0.05±0.03 


0.04±0.02 




SAAO A 


0.18±0.03 


0.17±0.03 


0.17±0.02 




SAAO B 


O.OOiO.ll 


0.00±0.11 


0.00±0.10 




SAAO C 


0.06±0.03 


0.06±0.03 


0.06±0.03 




OGLE 






0.06±0.01 


9v 


Perth 


0.14±0.03 


0.15±0.03 


0.15±0.02 




Canopus A 


0.47±0.13 


0.20±0.13 


0.20±0.13 




Yale 


0.12±0.08 


0.13±0.08 


0.13±0.07 




SAAO A 


0.30±0.05 


0.23±0.05 


0.23±0.04 




SAAO B 


0.00±0.33 


0.00±0.33 


0.00±0.33 




SAAO C 


0.00±0.07 


0.00±0.07 


0.00±0.07 


x 2 




729.7 


565.3 


720.9 


# points 




600 


595 


595+159=754 


xVd-o.f. 




1.28 


1.00 


1.00 



Table 2 

Point-source Point-lens Fit Parameters 



"PLANET "cleaned" data set; includes all data. 

^PLANET "super-cleaned" data set; does not include outliers with residuals > 3a. 
c Combincd PLANET "super-cleaned" and OGLE data sets. 



Mass Ratio 


AXthrcsh — 25 


A Xthrosh = 50 


Ax t \ rcsh = ioo 


10 u.u 


0.19-4] a 


0.19-4] 


0.19-4] 


1Q -0.5 


0.19-4] 


0.19-4] 


0.19-4] 


10 -i.o 


0.19-4] 


0.26-4] 


0.28-4] 


10 -1.5 


0.28-4] 


0.29-3.86 


0.29-3.48 


1Q -2.0 


0.38-2.87 


0.39-2.63 


0.40-2.37 


1Q -2.5 


0.50-1.96 


0.57-1.73 


0.67-1.51 


10 -3.0 


0.82-1.27 


0.85-1.20 


0.88-1.05 



Table 3 
Excluded Separations 



"Bracket indicates that companions are excluded at the largest separation we calculate. 



Mass Ratio 


A x " = 25 


A X " = 50 


Ax" = 100 


10 o.o 


1.00 


1.00 


1.00 


10 -0.5 


1.00 


1.00 


1.00 


10 -i.o 


1.00 


1.00 


1.00 


lQ-1.5 


1.00 


1.00 


1.00 


1Q -2.0 


1.00 


1.00 


1.00 


1Q -2.5 


1.00 


0.99 


0.98 


10 -3.0 


0.94 


0.87 


0.80 


1Q -3.5 


0.66 


0.54 


0.44 


1Q -4.0 


0.33 


0.24 


0.19 


1Q -4.5 


0.14 


0.11 


0.08 


1Q -5.0 


0.05 


0.03 


0.02 



Table 4 

Lensing Zone Detection Efficiencies 



